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Abstract 

We study a holographic gauge theory living in the AdS 4 space-time at finite tempera¬ 
ture. The gravity dual is obtained as a solution of the type IIB superstring theory with 
two free parameters, which correspond to four dimensional (4D) cosmological constant 
(A) and the dark radiation (C) respectively. The theory studied here is in confining 
and chiral symmetry broken phase for A < 0 and small C. When C is increased, the 
transition to the deconfinement phase has been observed at a finite value of C/\\\. It 
is shown here that the chiral symmetry is still broken for a finite range of C/ |A| in 
the deconfinement phase. In other words, the chiral phase transition occurs at a larger 
value of U/|A| than that of the deconfinement transition. So there is a parameter range 
of a new deconfinement phase with broken chiral symmetry. In order to study the prop¬ 
erties of this phase, we performed a holographic analysis for the meson mass-spectrum 
and other quantities in terms of the probe D7 brane. The results are compared with a 
linear sigma model. Furthermore, the entanglement entropy is examined to search for 
a sign of the chiral phase trantion. Several comments are given for these analyses. 
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1 Introduction 


The holographic approach [HEIIS] is expected to be applicable also to the supersym¬ 
metric Yang Mills (SYM) theory in curved space-time as well as in the flat Minkowski 
space-time. It would be interesting to make clear the properties of SYM theory in the 
curved space-time and in cosmologically developing universe. In this direction, some 
approaches have been extended to the field theory in the Friedmann-Robertson-Walker 
(FRW) type space-time [4J, [5] [6] . In this case, the 4D cosmological constant (A) can 
be introduced as a free parameter in obtaining the 5D sector of the 10D supergrav¬ 
ity solution. The dynamical properties of the 4D SYM theory on the boundary are 
then characterized by the sign of A. Through the holographic approach, it has been 
found that the SYM theory is in the confinement (deconfinement) phase for negative 
(positive) A |IJ|5J[6]. This implies that the dynamical properties of the SYM fields are 
largely influenced by the geometry of the background space-time. 

Furthermore, in this approach, one more free parameter ( C ) can be introduced as 
an integration constant in the solution of the supergravity. At first, this term has been 
added as the ’’dark radiation” to the 5D supergravity solution in the context of the 
brane world ®E\. Since it is defined in the 5D space-time, then its meaning in the 
4D theory was mysterious. In this context, afterward, it could be interpreted as the 
projection of the 5D Weyl term [9], TO]. On the other hand, from the viewpoint of 
holography, it has been found for A = 0 that C corresponds to the thermal radiation 
of the SYM fields at a finite temperature mi, then the system is in the deconfinement 
phase. Actually, the 5D metric in this case can be rewritten to the AdSs-Schwarzschild 
form, where C(> 0) corresponds to the black hole mass in this metric. 

It is easy to imagine that the above two parameters, C and the negative A, compete 
with each other to realize the opposite phase of the theory, namely the confinement 
and the deconfinement respectively. In fact, we find a phase transition at the point 
where these two opposite effects are balanced [IT, T2, JT. mm- A s a result, the SYM 
theory is in the deconfinement phase for b 0 > r 0 . Here the density of dark radiation 
C and the magnitude of | A | are denoted by using b 0 and r 0 , which are shown in the 
formula (12.5jl of the next section. By using these parameters, the phase diagram of the 
SYM theory in the FRW space-time is obtained as in the Fig. [I] 

In the deconfinement phase, the temperature T, which is given by the Hawking 
temperature of the 5D metric, appears. The critical temperature (T c ) of confinement 
deconfinement transition is given as T c = 0, which corresponds to the critical line 
r 0 = bo in the Fig. [T] In the region bo > ro, the temperature monotonically increases 
with bo for fixed r 0 . 

It becomes possible to study the properties of the finite temperature SYM theory in 
the deconfinement phase in the FRW space-time, where the three space is hyperbolic 
one. As mentioned above, for the case of A = 0 and finite bo, the bulk solution is 
reduced to the well known AdSs-Schwarzschild. So the SYM theory in its special limit 
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Fig. 1: The line (a) shows the critical line ro = &o separating the quark-confinement 
phase (A) and the deconfinement phase (B), (C). The lower critical line (b) ro = 0.76&o 
(%-br) is separating the chiral symmetry broken phase, (A), (B) and the symmetry 
restored phase (C). 


of this background has been studied already, see for example [16j. In this case, the 
two critical temperatures of confinement-deconfinement (T c ) and the chiral symmetry 
restoration (T x ) transitions are the same, namely T c = T x = 0. On the other hand, 
we show here that the value of T x shifts from T c = 0. Namely we find 0 = T c < T x in 
the case of finite A(< 0), and then the critical line is given by bo = 0.76ro. It is shown 
in the Fig. |T|in the parameter space of b 0 and r 0 . Therefore there exists a new phase 
in the region B of the Fig. CD where quark and gluons are not confined but the chiral 
symmetry is not yet restored. 

This implies a non-trivial thermal property of the SYM theory in the FRW space¬ 
time for negative A. In order to show and understand the details of this point, we 
have examined the embedding of the D7 brane in the background. For the embedded 
D7 brane, we have examined also the spectrum of the Nambu-Goldstone boson and a 
massive meson modes in the region of the newly found deconfinement phase where the 
chiral symmetry is still broken. Then we compare our results with some theoretical 
consequences obtaind from some typical phenomenological models for the chiral sym¬ 
metry breaking of QCD in order to make clear the dynamical properties implied by 
our holographic model. 

The outline of this paper is as follows. In the next section, a 5D space-time is given 
by a solution of supergravity as the dual of SYM theory in the AdS 4 background. 
In Sec.3, the spontaneous chiral symmetry breaking is studied by embedding the D7 
probe brane, and then the spectra of two scalar mesons are examined. Since one of 
them corresponds to the Nambu-Goldstone boson, then the Gell-Mann-Oakes-Renner 
relation and pion decay constant are also examined. In Sec.4, the mass relations of the 
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two scalars are discussed in terms of the sigma model. We suggest a modified sigma 
model which implied by our holographic analysis. In Sec.5, we searched a sign of the 
chiral transition in the entanglement entropy. We could find a change of the behavior of 
the entanglement entropy, but it could not lead to the phase transition. The summary 
and discussions are given in the final section. 


2 Gravity dual 


The holographic dual to the large N gauge theory embedded in 4D space-time with 
dark energy and ” dark radiation” is solved by the gravity on the following form of the 
metric, 

r 2 - Ft 2 

ds \ 0 = — (— n 2 dt 2 + A 2 al(t)'y i j(x)dx l dx j ) -I—- dr 2 + R 2 dFl\ , (2.1) 

where 

7 ij(x) = Sij 7 2 (x) , 7(x) = 1/^1 + > f2 = > ( 2 - 2 ) 

and k = ±1, or 0. The arbitrary scale parameter Tq of three space is set hereafter as f 0 = 
1. The solution is obtained from 10D supergravity of type IIB theory [TT, 12| 1131 14j. 

The factors A and n are obtained by introducing two free parameters as mentioned 
below. Here, we use the following form of solution, 


A = 


n = 



(2.3) 

(2.4) 


r 0 =v77?V2. b 0 = Rc 0 , co = Cfl 3 /(4oJ). (2.5) 

where the dark radiation C is introduced as an integral constant in solving the equation 
of motion. This solution expresses the case of negative A. Here, the ’’dark energy” (or 
cosmological term) A(t) is introduced by the following equation, 


^V + A = a 

\ a o/ a o 

Although it is possible to consider a time dependent A as in [m, we set it here as a 
constant A for simplicity. In the following, our discussion would be concentrated to the 
case of negative constant A and we assume small time derivative of ao(£). We should 
notice the following fact that the solution a 0 = 1/ y/jX[ =constant is actually allowed 
for negative constant A and k — — 1. 


( 2 . 6 ) 
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Chiral phase transition at finite temperature 


3.1 D7 brane embedding 

We study the chiral condensate and the q — q meson spectrum of the boundary theory 
by embedding the probe D7 brane for the flavor quarks. The D7-brane action is given 
by the Dirac-Born-Infeld (DBI) and the Chern-Simons (CS) terms as follows, 

S D 7 - -7- 7 I det (g ab + 2 Ta'F ab ) 

+T 7 /£(e 2 ' aV ™ A(, mI ) o , (3.1) 

i 

g ab = d a X< J d b X"G tw , c (Q ,.. ai) = d ai X^ ...d^C^ , 

where 77 is the brane tension. The DBI action involves the induced metric g ab and the 
U(l) world volume field strength F( 2 ) = dA^y 


The solution given in the previous section is obtained for the case of constant dilaton. 
So it does not play any role in the present case. For simplicity, we consider the dual 
theory on the boundary at r = 00. Then the induced metric for the above D7 brane is 
obtained as follows. Consider the above background (12.11) and rewrite it as follows 


ds 2 

us 10 


r 2 

—ds 2 


R 2 


( 4 ) 


dr 2 + R 2 dQ 2 


5 > 


r 2 , 2 R 2 

Jr- ds ^ + P? 


dp 2 + p 2 dni + dx ^ 


where 


(3.2) 

(3.3) 


ds 2 4 ) = (—n 2 dt 2 + A 2 al(t) , yij(x)dx l dx^ , (3.4) 

r 2 = p 2 + (X 8 ) 2 + (X 9 ) 2 . (3.5) 

Then, the induced metric of the D7 brane is obtained as 

ds l = Jp ds h + “T (( 1 + w ' 2 ) d P 2 + P 2 ^s) i (3-6) 

where the profile of the D7 brane is taken as (W 8 , X 9 ) = (w(p), 0) and w' = d p w, then 


2 2 1 2 

r = p + w . 

In the present case, there is no R-R hied, so the action is given as 

S D7 — — r 7 0 3 J d A xal(t)rf(x) J dpp 3 nA 3 ^Jl + w' 2 (p ), 


(3.7) 


(3.8) 
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where denotes the volume of S' 3 of the D7’s world volume. 


From this action, the equation of motion for w is obtained as 


w" 



The constant w ^ 0 is not the solution of this equation, so the supersymmetry is 
broken except for the case of trivial solution w = 0. 

3.2 Embedded solutions and chiral symmetry breaking 

In the confinement region, the chiral symmetry is spontaneously broken as shown in 
HU- The phase transition of the present model to the deconfinement occurs when 
the density of the dark radiation increases and excesses a critical point, which is 
given by 6 0 = r 0 (See the Fig. [[]). From this point, the Hawking temperature, 
Th(— ^r-^/l — (r 0 /&o) 2 ), appears for b 0 > r 0 . At this stage, when the tempera¬ 
ture appears, the chiral symmetry is usually restored. In the present case, however, we 
could observe the chiral symmetry restoration after transfering into the deconfinement 
phase at a finite which depends on r 0 , as shown in the Fig]]] above. These facts 
are explained below through the numerical analysis. 

All solutions of the above equation (13. 9 j) have the following asymptotic form 



(3.10) 


at large p. In the second term of the right hand side of (13.101) , the term proportional to 


log(p) arises from the loop corrections of the SYM theory since the conformal symmetry 
would be broken due to the existence of the cosmological constant in the present case 
HIEIE]. We could show however that this term is proportional to the quark mass m q , 
which is given by the asymptotic value of u’(oo). In order to see the spontaneous chiral 
symmetry breaking, it is enough to see chiral condensate (TT) = c = p 2 w\ p ^ t00 . Then 
the analysis is simply performed for m q = 0. 

We have therefore examined the numerical solutions for m q = 0 at various points 
in the parameter space bo — tq in order to find the transition points. Three typical 
solutions are shown in the left figure of Fig. [2] They are classified as Minkowski type 
(M-type) (u’(O) > r H = \Jb q — Tq), Black-hole type (B-type), which ends on r H at 
P = Pmin, and the trivial solution w = 0, which is always the solution of Eq. (l3.9j) . 

We notice that the above M and B types of embedded configuration of D7 brane 
correspond to the one of connected and disconnected D8 and D8 branes in the Sakai- 
Sugimoto model at finite temperature D2I 
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Fig. 2: Left; The solutions of w(p) with m q = 0 for bo = 1.3ro, and ro = 1.0. The 
solution (c) represents the trivial one, w — 0. The dashed curve denotes the horizon 
r H = 0-83). Right; The c — m q relations of embedded solutions for 

bo = 1.2,1.3, and 1.4 are shown with the same other parameters of the left figure. For 
the case of b 0 = 1.3, the embedded solutions at the three points, (a), (b), and (c), are 
shown in the left figure. 


We performed the numerical analyses for fixed ro = 1 by varying bo. In this case, 
the region of bo is separated to the following three one. 

• For b 0 > 1.31, we find only the trivial solution for m q = 0. Then the chiral 
symmetry is restored. 

• For 1.28 < b 0 < 1.31, there are three types of solution mentioned above. 

• For bo < 1.28, there are M-type and the trivial solutions. 

The distributions of the general solutions including these three types are shown in the 
right of the Fig. [2] for bo = 1.2, 1.3, and 1.4. Then one might consider that the chiral 
symmetry may be broken in the region b 0 < 1.31, however we must compare the free 
energies in order to see which solution is favored for the given b 0 when plural solutions 
for the same m q exist. 

The free energy for each solution is obtained by substituting the solution w(p) into 
the Wick rotated Eucledian D7 action (13.81) . Here the normalized free energy, 

E = j dpp 3 nAA\j 1 + w' 2 (p ), (3-11) 

is evaluated, and then the values subtracted the one for the trivial solution are shown 
in the Fig. [3] 

This figure shows the followings. 
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Fig. 3: Left; The blue (red) curve denotes the free energy difference, E — trivial, 
between the one for the Minkowski (Black Hole) embedding solutions of w(p) with 
m q = 0 and for the trivial solution w = 0. Right; The part, where the blue curve 
crosses zero, of the left hand figure is enlarged. 


• The free energy of B-type solution is always larger than the other two. Then this 
type of solution can not be realized. 

• The value of EM-type ~ -^trivial crosses zero at b 0 = 1.31. This implies that a 
transition from M-type solution ((TT) = c > 0) to the trivial solution ((TT) = 0) 
occurs at this point. This point is therefore the chiral phase transition point. The 
order parameter (TT) has a gap at this transition point. 

The transition line b 0 = 1.31 r 0 is shown in the Fig. [Q As a result, we find two critical 
lines (a) and (b) in the parameter plane (6o,r 0 ). The line (a) represents the transition 
point from confinement to deconfinement and (b) does represent the critical line from 
chiral symmetry breaking phase to the restoring phase. This result implies the fact that 
the density of the dark radiation necessary for the restoration of the chiral symmetry 
is larger than the one needed for realizing the deconfinement phase. 

We remember that the role of the dark radiation is to screen the force needed for 
the confinement. The same kind of force would be necessary for the spontaneous mass 
generation of massless quarks. The above result, that the chiral transition needs larger 
value of C than the case of the confinement- deconfinement transition, implies that the 
range of the force necessary to break the chiral symmetry is shorter than the one for 
the confinement. 
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4 Chiral phase transition 

and the Nambu-Goldstone (NG) bosons 

We study the meson spectra by solving the equations of motion for the fluctuation of 
fields on the D7 brane. In the A phase in Fig.l, we find the Nambu-Goldstone (NG) 
bosons and the Gell-Mann-Oakes-Renner (GOR) relation [18] for massless and small 
mass quarks, respectively. While in the B phase, in spite of the fact that quarks are 
deconhned, we expect the existence of the NG boson and a massive scalar mode due 
to spontaneous chiral symmetry breaking. 


To compute the meson spectra, let us consider the fluctuations of the scalar fields 
defined as 

x 8 = w(p) + 4> 8 , x 9 = 4> 9 . 

(p 8 and q i 9 correspond to the NG boson and a massive scalar boson, respectively. Note 
here that this phenomenon also represents the breaking of a global 17(1) symmetry. 
The wave functions are given in the following factorized form 

i h = f k (t ,i‘)^(p)y,(S 3 ), (4 = 8,9) 

where W(*S' 3 ) denotes the spherical harmonic function on three dimensional sphere with 
the angular momentum l. 


Then the linearized field equations of <f> 9 (p) and 4>f(p) for w ^ 0 are given as follows 


1 


dp&i + ~r^p{^o)d P (t>i + (1 + w' ) 


,R 4 mg 1(1 + 2) 


r n 2 


2 K, 


L 0 = p 3 nA 3 


\/l + 


w 


/ 2 ■ 


*<U = 


1 

nA 3 


(i) 


d r 2 (nA 3 


*? = 0, (4.1) 

(4.2) 


and 


9 p 4>i +1 -1 dp(Li)d p (l) l + (l+w/ ) 
-A 


- At7 - 2(1 + w' 2 )(A- (1) + 2'Ui 2 /\(2)) 


r n 




= -2 —d p (L 0 w w'K (l) )q if 




1 + w' 2 ’ 


Km = 


nA 3 


(9^2 (nA 3 ). 


(4.3) 

(4.4) 


Here are some remarks. 

i) For the 4D part of the wave-function, <p fc (a7), we have assumed the following 
eigenvalue equation, 

= = mlv k (x») 

V 94 


n A ip k (x») 


(4.5) 














where g 4 = —detand 


g^dx^dx" = (— dt 2 + al(t)'f i j(x)dx l dx j ) . 

(4.6) 

ii) The operator IZI4 is derived from the eight dimensional Laplacian □§ 
()3.6j) induced on the D7 brane. In fact, D 8 is expanded as 

for the metric 

□s - _ d a g ab Pm d b , 

Vm 


= + ■ ■ •, 

Vm 

(4.7) 

Vm = Vmp 6 P l + w ' 2 > y/m = PoiAyn, 

(4.8) 


where a, 6 = 0 ~ 7, p, ^ = 0 ~ 3 and the elliptics denotes the derivative terms with 
respect to the other coordinates, p and the one of S 3 . 

By using the above expansion and the approximation to neglect the time derivative 
of ao(t), we arrive at the above equations (14.ip and (14.31) . which are used to hnd the 
meson spectra (6j. 


4.1 Numerical Results 


Fig. 4 is the numerical results for the quark mass dependence of the mass eigenval¬ 
ues ms(m s ) and 7779 ( 777 ^ 0 ), which correspond to 0 8 (massive scalar) and 0 9 (Nambu- 
Goldstone boson), respectively. From the left panel of Fig. 4, we hnd the NG boson 
mass behaves as m 2 NG oc m q , which is consistent with the GOR relation. 

The GOR relation is expressed as 


777 


2 _ 

NG ~ 


2777 ? (TT) 
J 7T 


(4.9) 


where f n denotes the pion decay constant. Since (TT) is obtained for various 60 s 
through the solution of the D7 brane profile, then we can determine the 60 dependence 
of f n through the GOR relation (14.9)1 . Fig. [5] is the numerical results. 

Here is our findings. f n is almost constant within the confinement regime (60 < r 0 ), 
and then it increases with 60 in the deconfinement regime ( 6 0 > r 0 ). This result seems 
to be reasonable since the decay channels may increase in the deconfinement regime. 


4.2 Comparison with the linear a model 

In principle, it would be possible to derive an effective theory of mesons from the D7 
brane action as a functional of 0 8 and 0 9 . Then the correspondence of the parameters 
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Fig. 4: Left:The quark mass(m g ) dependence of m 2 NG (Nambu-Goldstone boson (</> 9 )). 
Right: The m 2 s { red line) for massive scalar(0 8 ) is shown by the solid lines. The dotted 
lines are the prediction from the sigma model given below. The line (1) denotes the 
one for (14.12j) and (2) for (I4.16I) . 
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Fig. 5: bo dependence of the chiral condensate (TT) and the pion decay constant f n 
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in the D7 brane action to the one of the sigma model will be obtained. However, there 
are various possibilities for higher order terms of the meson fields,. Then the resultant 
effective action would get very complicated, which is not useful to analyze. 

Instead let us compare the results we obtained in the above with those from the 
linear sigma model. At the mean held level, the Lagrangian density of the linear sigma 
model is given by 

£ = y(a 2 + vr 2 ) - ^((j 2 + vr 2 ) 2 + ha, (4.10) 

where /i, A and h are the parameters while a and n are the mean fields. The last 
term proportional to h, which plays a role of the quark mass term, breaks 17(1) chiral 
symmetry explicitly. For the small value of h, the vacuum is determined from the 
stationary conditions = 0, which lead us to 

(o-, tt) = (/*+2^, o) (4.11) 

/ 2 

with f n = y Then one obtains the mass spectra such that 

Ml(h) = A Ml(h) = 2fi 2 + A = 2^ 2 + 3 Ml(h) (4.12) 

J tv Jn 


In the right panel of Fig.4, we compared the results (I4.12|) with those obtained from 
the holographic method. Note here that the parameter region for bo is restricted to the 
confinement phase (6o < To)- Hereafter, we set ro — 1. The reason why we consider the 
linear sigma model only for the confinement phase is obvious because the Lagrangian 
density is solely written in terms of mesons. 

The results (I4.12p definitely deviates from our holographic ones. It rather looks to 
have a good fit to the result with b 0 = 1.3. So we might modify the sigma model in 
order to obtain a better fit to our results in the region of bo < 1. 

This is performed by adding the next order term of (a 2 + i r 2 ) as follows: 

c = k(a 2 + it 2 ) - + n 2 ) 2 - L(a 2 + a 2 ) 3 + ha. (4.13) 

2 4 o 

This type of modification is justified as far as it is based on the derivation of the 
effective sigma model by expanding the D7 brane action in terms of the n and a fields. 

The vacuum is determined as the configuration (a, it) = (a 0 , 0), which is the real 
solution of the following equation: 

h = (—/i 2 + Aiag + A2ag)ao . (4-14) 
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M 2 — —/i 2 + 3 Ai<7q + 5A 2 ctq • (4.15) 


Then the meson masses are given as 
Ml = -/i 2 + Aio-q + A 2 cto = —, 

From (14.151) . we find 

Ml = 2 Ml + p 2 

by setting the following relation among the parameters, 

Ai 


A 2 — 


3a 2 


(4.16) 


(4.17) 


This relation is a sort of tuning to obtain (I4.16j) . which provides a better fit to our 
holographic results (see the right panel of Fig. |4j. 


5 Entanglement Entropy near the transition region 


We study the entanglement entropy (See) near the transition region to find a signature 
of the phase transition. See is given by calculating the minimum area of the surface 
A whose boundary dA is set at the boundary of the bulk. As given in [ 231 [24] the 
holographic entanglement entropy is given by 


(~i _ S area 

= ^5j’ 


4 G 


N 


(5.1) 


where S area denotes the minimal surface whose boundary is defined by dA and the 
surface is extended into the bulk. 

We see the regularized finite part Sfinite- The detailed calculations are given in 
mm- This quantity contains two contributions from the curvature ro and the dark 
radiation b 0 . In order to see how the dark radiation affects on the entropy, we consider 
the following quantity, 

S finite = Sfi n ite Sfinite | bo=0 i (5 • 2) 

by subtracting the S finite \ bo=0 from S finite . 

Fig. [6] shows the relation between Sfinite/V and bo. V is the volume of the sphere 
with radius p = p 0 in FRW 4 space, 


V 


47raQ 



27raQ 


( 4po(4 + Po) 

V (pg-4) 2 


+ log 


2 ~Po \ 
2 + Po) 


(5.3) 


From the results shown in the Fig. [6l we can observe the followings. 


• As shown in the Fig. [61 for small b 0 region, Sf in it e /V is small and Sf in it e /V oc b^. 
The energy density of the dark radiation is also proportional to b$. Thus the 
small deviation of Sfinite and the energy density will lead to the first law of the 
thermodynamics. On this point, we will discuss more in the near future. 
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Fig. 6: The left figure is the relation between Sfinite./V and bo for r 0 = R = 1 and 
a 0 = 0.5. For large b 0 region and small b 0 region, Sf inite /V can be fitted by Sfi nite /V = 
26.1 Tjj (blue line) and Sfi nite /V = 0.976o(red line) respectively. The right figure shows 
the relation at small bo region. There is no specific phenomenon at the chiral transition 
point( bo = 1.31) 


• For large bo > 1 region, on the other hand, it increases with bo rapidly and 
Sfinite/V oc Tfj where the Hawking temperature Tf/(6 0 ) is given as [15] . 

r„(6„) = ^yi-(r 0 /6„) 2 . (5.4) 

Thus, at large bo region, the entanglement entropy becomes the thermal entropy. 
This behavior is expected as a high temperature behavior of the entanglement 
entropy S fin ite/V. 

• However, as far as we observe the Fig. O there is no specific signature at the 
chiral transition point (bo ~ 1.31) in the behavior of the entanglement entropy. 

6 Summary and Discussion 

We have studied SYM theory in the AdS 4 space-time. The holographic dual is expressed 
by a 10D supergravity solution which is described by two free parameters corresponding 
to the negative 4D cosmological constant and the dark radiation. These two quantities 
work to opposite direction to realize a typical phase of the theory. The negative A leads 
the theory to the confinement phase. However the dark radiation prevents it. Then we 
find confinement-deconfinement transition at their balanced point, r 0 = b 0 , as shown 
previously. 

Here we have pointed out that the chiral symmetry restoration does not occur yet 
at ro = bo as expected in the case of usual QCD. The chiral transition is found at 
bo = 1.31r 0 after the deconfinement transition. So there exists a new phase, where 
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chiral symmetry is broken but the quarks and gluons are deconhned. It is shown as 
the region B in the phase diagram in the Fig. [D In this region, the Nambu-Goldstone 
boson is certenly observed, and then we could examine for the mass spectra of mesons 
made of massive quark and anti-quarks to assure the GOR relation. We could study 
the mass relation of the NG boson and the massive scalar modes as expected from the 
sigma model, which describes well the spontaneous chiral symmetry breaking of QCD. 
While the modified sigma model might be consistent with our holographic results in 
the confinement region A, we could not find a simple sigma model which is consistent 
in both the regions A and B. 

Finally, we have examined the entanglement entropy to see the role of the dark 
radiation in the phase transition. In order to make clear the contribution of b 0 , the 
deviation of the entanglement entropy given at bo > 0 from the one at bo = 0 is 
numerically studied. For small &o(< r o)> the deviation is small and increases as 6 q, 
which is proportional to the energy density of the dark radiation. In the large 6 0 (>> r o) 
region, it increases rapidly and it is proportional to 6g oc T^. This behavior is the 
usual thermal behavior expected in the infrared limit of the theory. However a sharp 
signature of the phase transition has not been observed. 

We would like to mention about relations between our results and those from some 
4 dimensional effective theory. The Polyakov-Nambu-Jona-Lasinio (PNJL) model has 
been often used to obtain the QCD phase diagram at finite temperature/chemical 
potential. In this model, in addition to chiral symmetry breaking/restoration due to the 
quark-antiquark condensate which was originally developed in the Narnbu-Jona-Lasinio 
model, the confinement-deconfinement phase transition can be taken into account by 
the Polyakov loop potential. According to the model calculation, T x gets higher than 
T c . This is the similar result with ours. 

In the lattice QCD calculation, however, the result is the opposite. T x gets lower 
than T c [19] • In [20], the authors have tried to reproduce the same result from the PNJL 
model by including some extra terms, but they did not succeed. As the consequence, 
even among the 4 dimensional approaches, there is no common picture for the properties 
of QCD phase transition at finite temperature. Physically speaking, in the case with 
T x < T Cl massless baryons appear in the temperature regime T x < T < T c and they 
would affect the thermodynamic quantities such as the equation of state. 

On the other hand, the chiral phase transition in curved space has been recently 
discussed [23]. It is suggested that the negative curvature R < 0 (the negative cosmo¬ 
logical constant) shifts the critical temperature T x to the increased one. 
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